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Abstract 

A simple analytic tool namely the general series expansion method 
is proposed to find the solutions for nonlinear differential equations. By 
choosing a set of suitable basis functions {e n (t, to)}^= such that the so- 
lution to the equation can be expressed by u(t) = J^n^o c ne n (t, to). In 
general, to can control and adjust the convergence region of the series 
solution such that our method has the same effect as the homotopy anal- 
ysis method proposed by Liao, but our method is more simple and clear. 
As a result, we show that the secret parameter h in the homotopy anal- 
ysis methods can be explained by using our parameter to- Therefore, 
our method reveals a key secret in the homotopy analysis method. For 
the purpose of comparison with the homotopy analysis method, a typical 
example is studied in detail. 
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method; nonlinear differential equation 

PACS: 02.30.Mv, 02.60.Cb, 02.30.Hq, 02.30.Jr 



1 Introduction 

For a given nonlinear differential equation with initial and boundary conditions, 
it is well known that we can use the Taylor series expansion at the initial point 
to find its solution. On the other hand, according to the initial and boundary 
conditions, if we choose a set of suitable basis functions {e n (t, to)}n^o t° expand 
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the solution, we can obtain more efficient solution. For example, for a vibrate 
problem, we can choose periodic functions such as sine or cosine functions as the 
basis functions. Fourier series or transformation method is a typical example. 
In general, the convergence region of the series solution is relative small. How to 
enlarge the convergence region of the solution? The homotopy analysis method 
proposed by Liao[l] provides a way to solve the problem. We describe the Liao's 
method as follows. Consider the nonlinear differential equation 

Lu + N(u) = 0, (1) 

with some initial or boundary conditions, where L is a linear operator, and N 
is a nonlinear operator. The key step is to construct the following homotopy 

(l-p)[Lu- Lu ] + ph[Lu + N(u)} = 0. (2) 

When p changes from to 1, the solution of Eq.(2) changes from m to the 
solution of Eq.(l). Concretely, take the solution of Eq.(2) in the following series 
form 

+ oo 

u(t,p) = J2u n (t)p n , (3) 
n=0 

where u n (t) are unknown functions to be determined, for n = 0, 1, • • • . Sub- 
stituting u(t,p) into Eq.(2) and setting the coefficients of p n to zero yields the 
system of the differential equations of u n (t)'s. Solving the system gives every 
u n . In finally, taking p = 1 gives the solution u(t) = u(t,l) of Eq.(l). In a 
series of papcrs[2-10] and the book[ll], liao developed and applied his method 
to deal with a lot of nonlinear problems. According to different choice of uo(t), 
the final form of solution of Eq.(2) is given by 

m 

u(t)= lim y^Mm,n(/i)cne„(t), (4) 

m— >oo — ' 

n=0 

where {e n (t)}t=o * s a set °f functions and the approaching function fi m ,n{h) 
satisfies lmim^oo [i m . n (h) = 1. For example, if we take e„(i) = t n for n = 
0, 1, 2, • • • , we obtain the so-called generalized Taylor series solution. Liao points 
out that the generalized Taylor series provides a way to control and adjust the 
convergence region through an auxiliary parameter h such that the homotopy 
analysis method is particularly suitable for problems with strong nonlinear- 
ity[l,ll]. 

It is easy to see that the auxiliary parameter h is the key to the homotopy 
analysis method. However, the mathematical meaning of the parameter h is 
unknown so that the homotopy analysis method is still wrapped in a secret veil. 
This is an open problem in the theory of the homotopy analysis method. Thus 
a natural aim is to solve it. 

In the present paper, we propose a general series expansion method for non- 
linear differential equations. By comparing with the homotopy analysis method, 
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we give the answer of the problem on the mathematical meaning of h in the ho- 
motopy analysis method. In order to understand the key of my idea, we first 
prove that the so-called generalized Newton binomial series at the initial point is 
just the usual Newton binomial series expansion at another point to- For other 
basis functions, we also introduce an auxiliary point to whose mathematical 
meaning is clear. Indeed, we find that the point to can control and adjust the 
region of the corresponding series solution. Thus our method includes two key 
steps: one is the choice of the basis function, another is the choice of the point 
t - These two simple choices can provide us with an efficient way to expand 
the solution and control its convergence region. By the comparison of these two 
methods, we uncovers some secret aspects of the homotopy analysis method. Of 
course, as a simple but efficient method, the general series expansion method 
can be directly used to solve nonlinear equations. We will give other applications 
in the future. 



2 On the generalized Newton binomial theorem 

In order to understand easily the key of our idea, we discuss the Liao's gen- 
eralized Newton binomial theorem. For real number a(a ^ 0, 1, 2, 3, • • • ), the 
generalized Newton binomial theorem is given by 



(i+tr= iim 5>a ,n w t" 



(5) 



n=0 

whose convergence region is 



K*< y|-l,(-2</i<0), (6) 



where 



and 



/C n W = H0"- a £(-i) J '( „ )a + ^'. (7) 



n 



t(a — 1) • • • (a — n + 1) 



n nl 



(8) 



The approaching function n™' n (h) satisfies lim TO _ ) . cx) /i TOl „(/i) = 1 when n > 1. 
Liao[ll] points out that the generalized Newton binomial theorem provides a 
way to control and adjust the convergence region through an auxiliary parameter 
h. We next point out that the generalized Newton binomial theorem is exactly 
the usual Newton binomial expansion at another point to = — 1 — ^. Indeed, 
when t = —I — \, the Taylor expansion at the point to is given by 



n=0 
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= iim £/c*wu F < ^ 

where 

m — k / , \ 

WW = {-h) k - a r 7 ) (h + iy(-iy. (io) 

Correspondingly, the convergence region \t — to I < |1 + of the binomial 
expansion at the point to = — 1 — j- is just (6). This result reveals clearly the 
secret of the auxiliary parameter h and the essence of the generalized Newton 
binomial theorem. More details can be found in Ref.[T2]- Furthermore, we can 
prove that the generalized Taylor series is essentially the usual Taylor expansion 
at another point [IB]. 



3 General series expansion method 

For a given differential equation of function f(t) 

N(f,f, ••■) = 0, (11) 

with some original or boundary conditions at the point t = 0, e.g., /(0) = 
0>/'(0) = 0, /(+oo) = 0. Firstly, we choose a set of base functions eo(t, to), ei(t,to), • 
Then we expand the solution as a series , 

f{t) = ^a n e n [t,t )- (12) 

n=0 

Substituting the series solution into equation (11) and using the original and 
boundary conditions, we can determine the values of the parameters a n for 
n = 0,l,---. 

On the other hand, we can expand every e n (t, to) as 



m=0 

where 



e n (t,t ) = b nm (t )e m (t), (n = 0, 1, • • •), (13) 
e m (t) =e m (i,0). (14) 



Then we have 

/(f) = ^^] a n b nm (t )e m (t), (15) 

n— m— 

Therefore there is a relationship between the point to and the auxiliary param- 
eter h such that 

m 

/(f) = lim V//"<"(/i)c„ eil (i), (16) 

m— >+oo * — » 
n=l 
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where c n satisfies 

/(t) = ^c n e n (t). (17) 

71=1 

Remark 1. Comparing with the usual series expansion methods such as 
power series method whose expanding point is in general the original point, an 
important different point of our proposed general series method is that the ex- 
panding point of our series is a motive point to which can control the convergence 
region. 

Remark 2. Ma and Fuchssteiner proposed a powerful approach for finding 
exact solutions to nonlinear differential equations [14 J. The crucial idea is to 
expand solutions of given differential equations as functions of solutions of solv- 
able differential equations, in particular, polynomial and rational functions. A 
more systematical theory on decompositions and transformations is presented 
very recently in Refs.[TS] and [TB]- The resulting theory unifies many existing 
approaches to exact solutions such as the tanh- function methods, the homoge- 
neous balance method, the exp-function method and the Jacobi elliptic function 
method. 

Remark 3. The introduction of the parameter h in the homotopy analysis 
method is to control the convergence region. In order to increase the speed 
of convergence of the series solution, Liao used the Pade approximation to do 
it (see Ref.[ll] for details). Therefore, we can also use Pade approximation to 
increase the speed of convergence of the general series solution. 

4 A typical example 

In the chapter 2 of the book [11], Liao studies a typical example which reads 

V'(t) = l-V 2 (t), 7(0) =0, (18) 
whose exact solution is given by 

V(t) = tanh(f). (19) 
Liao constructs the following homotopy 

a- ^'^-^ l-^g^ + fft.)-!) (20) 
where q G [0, 1]. Then we have 0) = V (t), 1) = V(t). We take 

+ 00 

$(t;q) = V (t)+J2 V ™(t)q m , (21) 

m— 1 

where V m (t) are the functions to be determined. Substituting the solution into 
the homotopy equation yields a polynomial of q. Setting the coefficients of the 
polynomial gives a series of linear differential equations which solution V m (t) 
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can be obtained. In final, we take q = 1 to give the solution needed. According 
to the rule of solution expression, by taking Vo(t) — t, Liao obtains a generalized 
Taylor series solution 



V(t)= hm y2f^' n (h)a 2n+1 t 2n+1 , (22) 

where 



m — J-oo ■ 

n=Q 



"' n 'n-l + k 



k=0 



^> n (h) = (-hrj2 (i+Af, (23) 



and «2n+i are the coefficients of the Taylor expansion of tanh(i) at the point 
t = 0, 



tanh(t) = J2 a 2n+it 2n+ \ (24) 



n=0 



whose convergence radius is po ~ § (indeed, we have a 2 n+i = - — ( 2 fc+2)! — -^2fc+2 
where B 2 k+2 are Bernoulli numbers). Liao gives the convergence region of the 
generalized Taylor series (22) 



o<t< Po y / I |-i, (25) 

with —2 < h < 0. In the paper[T3], we have studied this case in detail and 
showed that the generalized Taylor series was just the usual Taylor series at 
another point. In the present paper, I study another set of basis functions. As 
the same as Liao [11], take the set of basis functions 

e n (t) = e - nt ,(n = 0,1,2, •••)■ (26) 

Liao gives the m-th order approximate solution 

V(t) « 1 + 2 exp(-2nt)/C»} - exp(-t){l + \ + £ exp(-2i)}™. 

n=l 

(27) 

When h = — 2, the solution becomes 

m 

V(t) « 1 + 2 ^{(-1)" exp(-2nt)} + (-l) m+1 exp(-(2m + 1)), (28) 

n=l 

which converges to the exact solution in the region < t < +oo. However, if 
we expand directly the exact solution using the basis functions {e n (i)}o*°°, we 
have 

V(t) = 1 + 2 X){(-l) n exp(-2n*)}, (29) 

n=l 



G 



which converges to the exact solution in the region < t < +00, and dose not 
converge at the point t = 0. This is a serious shortcoming for solution (29). Al- 
though the Liao's solution (27) can avoid this weakness, its accurate convergence 
region can't be given. Our method can overcome these two difficulties. 

We next solve the equation (18) using our general series expansion method 
at the point to to give the same results with Liao's. We take the basis functions 
e n (t,to) = {exp(— 2t) — t }™ an d expand the solution as 

+00 

V(t) = a„{cxp(-2i) - t }", (30) 

n=l 

and substitute it into Eq.(18) to give 

a = i^,a n = 2(-i rjYT L^,n = l,2,.... (31) 

Therefore the series solution is given by 

1 - t + °° 1 

n— 1 

+ 00 +00 ,k-n 

= l + 2VJ {E Q^ Tr} (- ir exp(-2„ t ) 

n— 1 k—n x / x ' 

k\ t n 



m—too 

n—1 k—n 



whose convergence region is 

- i ln(2t + 1) < t < +00, (33) 

which can be obtained by computing | a ^ +1 (exp(— 2t) — to) < 1 from (30). When 
to > such as to = 1, the convergence region consists of the region < t < +00. 
Thus t can control and adjust the convergence region. This shows that our 
method is more simple and efficient in some degrees than the homotopy analysis 
method. 

Take f = — 1 — ^ , and denote 

""•'"<") - 1 CO (dv* - <-"» n+i t 0+ ^ < 34 > 

k=n y ' k=n 

Therefore the solution can be written as 



V(t) = l+ lim 2V ^ m '"(ti)(-l)"exp(-2nt), (35) 

n=l 

whose convergence region is 

1 2 

- - ln(-l - -) < t < +00. (36) 
Zi n 

This result once again shows that the role of our parameter t is the same as 
the Liao's parameter h. 
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5 Conclusion 



We propose a general series expansion method to find the solutions to nonlinear 
differential equations. Because the expanding point of the corresponding series 
solution is a motive point, the method has a freedom to control the convergence 
region. Through detailed analysis of a typical example, we show that we can use 
the general series expansion to give the same result obtained by the homotopy 
analysis method. As a result, we reveal the secret of the auxiliary parameter h 
which is the key of the homotopy analysis method. In some degrees our method 
is more simple and efficient than the homotopy analysis method. 
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